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Abstract. - Many metastable complex fluids such as colloidal glasses and gels show distinct 
nonlinear viscoelasticity with increasing oscillatory-strain amplitude; the storage modulus de- 
creases monotonically as the strain amplitude increases whereas the loss modulus has a distinct 
peak before it decreases at larger strains. We present a qualitative argument to explain this 
ubiquitous behavior and use mode coupling theory (MCT) to confirm it. We compare theo- 
retical predictions to the measured nonlinear viscoelasticity in a dense hard sphere colloidal 
suspensions; reasonable agreement is obtained. The argument given here can be used to obtain 
new information about linear viscoelasticity of metastable complex fluids from nonlinear strain 
measurements. 



Rheological measurements are commonly used to investigate dynamical and structural 
properties of complex fluids. Oscillatory linear viscoelastic experiments measure the response 
function, or the complex modulus, G*(ui) = G'{lo) + iG"(u>) [1]. This provides information 
about the relaxation of the microstructure of the system as a function of oscillation frequency, 
lo. This function is strictly defined only within the linear response regime, and thus an 
important practical experimental issue is the determination of its range of validity. However, 
many complex fluids exhibit a nonlinear response even at very small applied strains. This can 
provide an alternate means of probing the properties of the system; the strain amplitude, I7I, 
can be used as an independent parameter [2,3]. While the nonlinear rheological properties 
of complex fluids are very diverse [1] , most colloidal suspensions and polymer solutions show 
a simple yielding behavior, where both the storage (G") and loss (G") modulus decrease 
monotonically as I7I increases. Less common is strain hardening (or strain-stiffening), where 
G' increases with |7|. This can be observed in polymer solutions with strongly cross- linked 
chains [1], soft-biological tissues [4,5], and in dense colloidal suspensions [6,7]. There is 
another type of observed rheological behavior that lies between these two kinds of responses: 
as I7I increases, G" continuously decreases above a critical strain, whereas G" exhibits a 
distinct peak before decreasing at even larger applied strains. This behavior is ubiquitous 
in metastable complex fluids. It was first reported in carbon black solutions [8]; many other 
systems exhibit similar behavior, including polymer systems [9], fumed silica [10,11], colloidal 
glasses [12,13], emulsions [14,15], gels [16] and electrorheological fluids [17,18]. Despite the 
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diversity of these systems, the observed behavior is strikingly universal: (i) The peak in G" 
appears around 1CP 2 < J — y ] < 1 over a wide range of frequencies; (ii) the peak in G" is 
observed for a frequency regime where the system is elastic, with G" > G"\ and (iii) the shear 
thinning exponent v of G" oc \'y\~ u at larger strains is always about one half of the exponent 
associated with G' . There have been many attempts to explain these universally observed 
phenomena. It has been argued that the increase of G" as a function of I7I is due to the 
increase of the effective volume of the temporal network structure in polymer solutions [9]. 
Alternative suggestions include the breakdown of aggregates into small clusters which are 
more dissipative [11] and the rearrangement of clusters due to the strain associated with 
thickening [17]. The phenomenological soft glassy rheology (SGR) model also accounts for 
the peak in G" as a function of I7I, and posits that flow near the glass transition arises from 
activated barrier crossing [19]. However, a full microscopic account of this ubiquitous behavior 
is still lacking. 

In this Letter, we show that the origin of the ^-dependent behavior of G* is in fact 
remarkably general and transcends the physical ingredients of the SGR picture; a peak in G" 
should be present in any viscoelastic complex fluid that behaves like an elastic solid at small 
strain but yields and flows at larger strains. In particular, we explicitly demonstrate that 
this is true even in the absence of local activated processes. To make this argument more 
quantitative, we focus on colloidal suspensions in the supercooled state and calculate the u 
and I7I dependence of G* using mode-coupling theory (MCT), generalized to account for the 
presence of oscillatory shear strain. This method successfully accounts for the peak observed 
in our experiments on a hard sphere colloidal suspension near its glass transition density. 

The presence of the peak can be qualitatively explained using the following simple argu- 
ment. Throughout the paper, an oscillatory strain, j(t) ~ \ j\e luJt , is considered. In the small 
strain limit, the stress is written as a(t) = G*(u)"/(t). As the simplest example, we consider 
a Maxwell model, where the storage and loss moduli are given by 

respectively. Here To is a characteristic relaxation time which in many systems represents the 
structural relaxation time. According to this equation, G"(lo) has a peak around lj ~ l/Vo- By 
contrast, for large strains many complex fluids experience a drastic decrease of their structural 
relaxation time. Such a decrease can be characterized by the strain rate amplitude, 7), and 
is given by the phenomenological expression 

1/t ^ 1/ro + K\i\ v , (2) 

where v is a positive exponent and if is a constant. If to is replaced by t in Eq. a 
peak of G" should be observed at to ~ 1/t ~ K\^\ v ] this equality is satisfied if u> ^> 1/to 
and the system's intrinsic relaxation is very slow (thus 1/to can be neglected). In particular, 
metastable systems which show slow dynamics, such as supercooled fluids and gels, have an 
exponent v very close to unity and K of order 1 [20,21]; then, the peak always appears 
near lo ~ 7I or I7I ~ 1. This argument also accounts for the facts that (i) the position 
of the peak is relatively insensitive to the frequency, (ii) the peak of G" should appear at 
frequencies where the linear viscoelasticity is dominantly elastic, (G' > G"), and (iii) the shear 
thinning exponents at larger strains are related and are given by G' ~ |7| " and G" ~ 7| _I/ . 
Considering the broad diversity of systems, it is remarkable that this simple argument correctly 
captures the features of the peak in G" . The above argument is based on the assumption that 
nonlinear effects solely stem from the strain amplitude I7I (or more precisely from the strain 
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rate amplitude |7| = Lo\"f\) and are independent of the time-dependence of the strain. In other 
words, the stress is given by a(t) — G*(lu, \"f\)"f(t). This is consistent with recent experiments 
which have taken the effect of higher order harmonics into account and shown that they 
account for at most 10% of the response [3,22], thus, lending support to the notion these 
terms are negligible. 

In order to make this argument more quantitative, we use a microscopic description based 
on MCT. First, we focus on a typical complex fluid which shows a peak in G" as a function 
of I7I, namely, a colloidal suspension in the supercooled state. To date, MCT is the only 
first principles theory that describes the dynamics of the supercooled state. Recently, several 
attempts have been made to extend the MCT approach to systems under shear strain for either 
stationary systems [20,21,23] or for transient states [24,25]. In all previous analyses, 7 has been 
assumed to be time- independent. Here we generalize the analysis [20, 21] to time-dependent 
oscillatory strain. We consider strain in the x direction given by "/(t) = \j\ sm(ujt) . The 
arguments used are exactly the same as for the constant shear case [20,21], which should be 
valid as long as the shear rate amplitude w|7| is small compared to the inverse of the diffusion 
time a 2 / D. We define a density correlation function which satisfies generalized translational 
invariance, F(k,t) = iV _1 (5/9 k (_ t )(f)(5p_k(0)}, where N is the total particle number and k(i) 
is the time dependent wavevector defined as k(i) = (k x , k y + j(t)k x , k z ) [20]. Then, the MCT 
equation for F(k, t) is 



where D is the diffusion constant in the dilute limit, k = |k|, and S(k) is the (unperturbed) 
static structure factor. M(k, t) is the MCT memory kernel and is a nonlinear function of 
F(k, t) and S(k), whose explicit expression is given in Refs. [20,21]. Solving Eq. numerically 
is a demanding task because of the anisotropic nature of the equation in the presence of the 
unidirectional strain. In order to make the computation tractable, we adopt the "isotropically 
sheared hard sphere model" (ISHSM) approximation proposed by Fuchs et al. [24] , where the 
anisotropy of the strain which appears in k(i) is neglected and the scalar |k(i)| is approximated 



by |k(i)| = ^Jk 2 + k x k y j(t) + k 2 j 2 (t) ~ ^/l + 7 2 (t)/3 k. This approximation is expected to 



be qualitatively correct as verified by the observation that the dynamics of supercooled liquids 
under shear strain is surprisingly isotropic [21]. We solve Eq. (|3J) using this approximation 
and evaluate F(k,t) for a given set of uj and |7|. The complex modulus is calculated in the 
same manner as M(k,t) in eq.© and is given by [21] 



where V(k) = k x S~ 2 (k)nS' (k) and n is the particle number density. We calculate G^ l (uj;j) 
for a hard sphere colloidal suspension with particle radius a for a volume fraction <j) = 4n7ra 3 /3 
near the MCT glass transition point, <fr c w 0.516, at which point F{k, t) undergoes an ergodic- 
noncrgodic transition [26]. In Fig. 1, we show the frequency and strain dependence of 
GJf L (oj, I7I). The volume fraction is set to <p — 4>c X (1 — 10~ 3 ). We calculate S(k) from 
the Percus-Yevick equation. The frequency is scaled by Aa 2 /D and G^ h is plotted in units 
of fcBT/8a 3 . The behavior of G^ L demonstrates that, in the small strain limit, a broad dip 
(f3- regime) is observed between 10~ 3 < lu < 10 1 . A peak at low frequencies u> ~ 10 -4 is the 
a-relaxation regime below which liquid-like relaxation is observed. In the /3-regime, G^ L is 
larger than G^ L and the system is solid-like. As the strain is increased in this regime, a weak 
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Fig. 1 - Frequency and strain dependence of (a) the storage modulus G^l^i 7) and (b) the loss 
modulus G^ L (w;7) for hard sphere suspension near the glass transition density. 



peak in G^ L is observed around | — ^ | ~ 0.4, whereas G' NL starts to decrease monotonically. At 
even larger strains, G^ L also starts to decrease. The decrease of Gj^ L is well described by a 
power law, G' NL ~ |7|~ 2iy and G^ L ~ |7| _!y with v = 0.9. The features that appear at both 
high frequencies and strain amplitudes occur in the regime where uj\^\ (or the Peclet number 
w|7|a 2 /D) becomes very large and thus the validity of the present MCT framework is ques- 
tionable. Thus, they should not be considered as physical. The position of the peak appears 
to be almost frequency independent. Although v varies with the system, these general features 
are shared by most experiments where the peak is observed [8-18]. It should be emphasized 
however that the peak in G^ L as a function of |7| is not generic, intrinsic shear thickening [5,7] 
but originates instead from a strong decrease of the structural relaxation (a-relaxation) time 
at large strains. Indeed, the a-relaxation time, t q , obtained from F(k,t) is found to behave 
as T a ~ |7| _!/ - Furthermore, our results show that the origin of the peak in G^ L is more 
general than that posited by the SGR model [19]. In particular, MCT does not capture the 
local rearrangements of particles that occur via activated barrier crossing, but it nevertheless 
properly captures the features of G^ L . We thus conclude that activated processes are not 
necessary to describe these features. 

While these calculations have been done for oscillatory strain, we note that the time- 
dependence of the generated strain, ^(t), has only a small effect on the ui and I7I dependence 
of GJj L . To show this, we calculate G^ L by approximating the oscillatory strain rate by a 
constant, ^(t) ~ u\"f\. GJj L thus obtained is qualitatively similar to that calculated for the 
case of oscillatory strain with a slight shift to a lower strain amplitude (I7I ~ 0.3) of the 
peak position. These conclusions support our qualitative argument based on the Maxwell 
model. The qualitative argument also asserts that G*(u>) is a function of Cj = lot (see Eq. flj) 
and the effect of the large strain comes in only through the relaxation time in the form of 
Eq. J2J ■ We find that this argument is valid for structural relaxation in glassy systems within 
the approximate framework of the MCT described in this paper. T a (ff) extracted from the 
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peaks of G^l shown in Fig.l (b) as a function of the strain rate is well described by Eq. J2J 
with To = t q (7 — > 0), v — 0.9 and a positive constant K (~ 2.5 in units of 4a 2 /D). We also 
find that, around the peaks, the shape of G^ L (a>;7) as a function of I7I matches well with 
G*(ujT a ) calculated at I7I = as a function of ur a , where the strain dependence of r a is given 
by Eq. Q. 

To demonstrate the power of the MCT approach, we make a quantitative comparison of 
theory with an experimental system, namely, a hard sphere colloidal suspension. Experiments 
are performed with sterically stabilized PMMA particles suspended in a mixture of decalinc 
and cycloheptyle bromide. Their radius is a=197 nm and their volume fraction is <f> « 0.56, 
which is a super-cooled fluid, just below the glass transition. Rheological experiments are 
performed in this supercooled state at room temperature using a strain-controlled rheometer 
(ARES, TA Instruments). Results for strain-sweep measurements at a fixed frequency w = 
0.2 rad/sec are shown in Fig. 2. A distinct peak is observed in G^ L at I7I « 0.15, while 
the onset of a monotonic decrease of G' NL is observed at the same strain amplitude. The 
shear thinning behavior is well described by a power law, with a shear thinning exponent 
v « 0.7. In the linear response regime where the strain is small, MCT is known to be 
successful at predicting the frequency dependence quantitatively and the overall shape of 
the complex modulus qualitatively near the glass transition points. Due to the generically 
approximate nature of MCT, a direct comparison to experimental data still requires several 
fitting parameters. It is known that MCT underestimates the glass transition volume fraction 
by more than 10% [27] (experiments predict the glass transition density around 4> g ~ 0.58, 
whereas MCT predicts the nonergodic transition at (f> c w 0.516). MCT also underestimates 
the amplitude of the complex modulus over a wide range of frequencies, as theoretically 
shown for a semi-dilute suspension [28,29]. Moreover, MCT entirely neglects the effects of 
hydrodynamic interactions. In order to determine the parameters to correct for the above- 
mentioned problems, we first perform a linear viscoelastic measurement for wide range of 
frequencies (10~ 2 < u> < 10 2 rad/sec) in the small-strain limit (I7I = 0.005 and 0.01). Then, 
we regard the volume fraction as a free parameter and determine it in such a way that the 
frequency dependence of GJj L (w, I7I = 0) calculated by MCT in the absence of strain matches 
with the shape of the measured G*(lu, I7 <C 1). The volume fraction thus determined is 
4> = 4> c x (1 — 10~ 2 ' 3 ). The hydrodynamic interactions are approximately taken into account 
by replacing the diffusion constant D in eq.© by Dm(k), the short time collective diffusion 
coefficient with full hydrodynamic interactions, which was calculated using the Beenakker- 
Mazur theory [30]. Finally, the amplitude of G* is matched with experimental data in the 
7 — > limit. For these parameters, we calculate Gj^ L for finite strains. The results thus 
obtained reproduce the peak in G^ L but overestimate the position of the peak. The lines in 
Fig. 2 are obtained by using the MCT calculations with a strain that is increased by a factor 
of approximately 3, so that the peak overlaps with the experimental data. The discrepancy 
of the peak position might stem largely from the ISHSM approximation which we introduced 
for technical reasons. Indeed, it is obvious that the ISHSM approximation underestimates the 
onset of shear thinning because it neglects the term linear in |7|. The shear thinning after 
the peak exponent from MCT is found to be v ~ 0.9 for G^ L and 2v for G^ L . While MCT 
overestimates of the shear thinning exponent and the height and position of the peak which 
is somewhat larger than that extracted from experiment, the overall shape and qualitative 
behavior of GJ^l is captured well. Considering the known problems associated with MCT even 
in the linear regime, the result shown in Fig. 2 generated by our "nonequilibrium" MCT is 
rather remarkable. 

Although in this Letter we have investigated a limited class of systems where the MCT 
approach may be applicable, we believe that our analysis is more general and gives significant 
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Fig. 2 - Comparison of MCT with experiment for a dense hard sphere suspension for uj = 0.2 
rad/sec. The empty and filled circles represent the storage modulus, G^li an d loss modulus, Gj^li 
respectively. Bold lines represents the MCT results using the ISHSM approximation near the glass 
transition density. The lines are shifted in the direction of the strain by a factor of 3. The modulus 
is plotted in units of k^T/8a . 

insight into the underlying physics of the peak in the loss modulus as a function of strain, 
observed in various types of complex fluids. The major origin of the peak in the loss modulus 
is the strong decrease of the structural relaxation time at large strains, and has no connection 
whatsoever with intrinsic shear thickening. This peak should always be observed for systems 
which exhibit slow dynamics and possess a broad band of solid-like behavior above a crossover 
frequency where the storage modulus is larger than the loss modulus. The theoretical argu- 
ment predicts that, near the peak of G" , the strain and frequency dependent viscoelasticity 
is well approximated by the linear viscoelasticity G* (ujT a (\j\)) , as a function of the frequency 
scaled by the relaxation time r a , whose strain dependence is given by Eq .(J2J. This implies 
that nonlinear measurements are a promising alternative for probing the structure and dy- 
namics of systems where linear viscoelastic experiments are limited to a narrow frequency 
window [31]. This will provide important insight into the physical origins of the nonlinear 
viscoelastic response of a wide range of soft materials. 

The authors acknowledge support from the NSF (# 0134969, KM and DRR), (DMR- 
0602684, HMW and DAW), the Harvard MRSEC (DMR-0213805) and the Swiss National 
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